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ABSTRACT 

Utilizing the Lorentz invariance of the correlation function we study the effects of 
anisotropic transverse flow on the HBT correlation function. In particular we show that 
directed flow would evidence itself by non zero "side-long" and "side-out" terms in the 
correlation function. We also show that the study of ratios of the correlation functions 
evaluated on different event subsamples and/or different Lorentz systems provides infor- 
mation on the source dimensions and velocity which could be less biased by corrections 
of different kinds. We present a fitting technique appropriate for HBT analysis which 
does not require binning and is especially useful for multidimensional fitting. 
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1 Introduction 



The Hanbury-Brown-Twiss (HBT) analysis of multiparticle production processes is 
becoming a widely used technique. It provides information on the space-time evolution 
of an excited strongly interacting system produced in high energy collisions. In recent 
years many experiments have collected high statistics data which permits detailed 
multidimensional analysis of the correlation function. From the theoretical point of 
view many valuable results have been obtained, both in the understanding the of the 
HBT method itself and in the ways it can be applied to the data 0—0. It is now 
clear, for example, that the correlation function depends sensitively on the dynamics 
of the expansion of the system. Most of the theoretical observations have been made 
using different kinds of models to describe this expansion. 

The discovery of transverse directed flow in the collisions of ultrarelativistic nu- 
clei |J implies that the HBT analysis of such data should take the effects of flow into 
account. In our previous publication [[| we used RQMD generated events to show how 
directed flow affects the extracted HBT radii. Utilization of the Lorentz invariance 
of the correlation function, defined as the ratio of the invariant two particle distribu- 
tion to the product of two invariant one particle distributions, can provide valuable 
information about the source and can be used as a check of the validity of the assump- 
tions contained in certain models. These questions are addressed in the current paper 
mostly in the context of directed flow. In particular we argue that directed flow would 
evidence itself by non zero "side-long" and "side-out" cross terms in the correlation 
function, even if the correlation function is determined with an azimuthally symmetric 
event sample (case in which the reaction plane is not determined). 

In dealing with the experimental data we find that the conventional method of 
fitting for the correlation function (especially multidimensional fitting) has several dis- 
advantages and sometimes gives biased results, in spite of many improvements and 
suggestions made recently The main problems relate to the low statistics in 

some bins and to a proper evaluation of fluctuations in the distributions of mixed 
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pairs. Multidimensional fitting is essential for the questions under consideration, and 
we present a fitting method without binning which is free from the problems mentioned 
above and is very flexible in its applications. 

The paper is organized in a following way. We start with Lorentz properties of 
the correlation function and discuss how the extracted source radii depend on the 
rapidity of the frame in which the analysis is done, and how anisotropic transverse flow 
affects the correlation function. We then discuss the ratio of the correlation functions 
evaluated with different event subsamples and/or in different Lorentz systems and 
discuss how these ratios could help in the selection of the appropriate functional form 
for the correlation function and in the extraction of parameters of the source such as the 
anisotropic flow velocity. Finally we present the new method of fitting the correlation 
function which does not require binning. 

2 The HBT correlation function in different frames 

The objective of an HBT study of the space-time evolution of multiparticle production 
is typically the determination of the effective source dimensions, time duration of parti- 
cle emission and the velocity of the source. Boson interferometry is a very effective tool 
for this purpose, but like any other tool has its own limitations. The first limitation, a 
rather trivial one, is that one can study only that part of the source which emits parti- 
cles into the experimental acceptance. This very simple statement mathematically can 
be written as space-momentum correlation of the source function. It implies that the 
study of the correlations between particles produced in different parts of phase space 
permits one to measure the source size as it seen from different directions. 

The second limitation is not so trivial. It arises from the physics of interferometry. 
Interferometry measures the distance between pions at the instant the second pion is 
produced. This distance depends not only on the size of the source, but also the duration 
of the emission (which affects the distance traversed by the first pion before the second 
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pion is produced) and the velocity of the source (the source could move to another 
place before the production of the second pion). Mathematically it is related to the 
fact that only three of the four components of the pair momentum difference g« are 
independent, because of the particles being on the mass shell. The fourth component 
of the momentum difference can be expressed through first three components and the 
pair velocity V: 

* = Vq. (1) 

Often this results in the situation when there are more unknowns than the number 
of parameters possible to extract experimentally, and additional assumptions must be 
made about the source to resolve the problem. 

The correlation function, being a ratio of invariant distributions, is invariant with 
respect to a Lorentz transformation: 

_ d 6 n/d 3 Pl d 3 p 2 _ 
6 (q ' P) " dV^PidVd 3 ** " (q ' } ' { ) 

where P = pi + P2 is the total momentum of the pion pair, and the prime denotes the 

values in any other Lorentz system. Lorentz invariance provides us with the possibility 
to calculate the correlation function in any frame provided we know its functional form 
and value in any other frame. Below we use the following notation: by an asterisk (*) 
we denote the values in the source rest frame (but we use y* for the source rapidity in 
the laboratory frame), and a tilde (~) is used for the LCMS (Longitudinally Co-Moving 
System) frame, in which the longitudinal velocity of the pair is equal to zero (V z = 0). 

For example, using this notation, the Lorentz transformations between the analysis 
frame (with rapidity y) and the source rest frames reads: 

q* z = q z cosh(y - y*) + q sinh(y - y*) , (3) 

<?o = Qo cosh(y - y*) + q z sinh(y - y*). (4) 

The longitudinal and transverse velocity components of the pair in the analysis frame 
are: 

V z = tanh(y -y), (5) 



V ± = V ± /cosh{y-y). (6) 
Very often the correlation function is written in the form: 

C(q) = 1 + Xe- Q{q) . (7) 

Using the small q approximation in the source rest frame (we discuss below the existence 
of such a frame ) one can represent Q as: 

Q*2 i *2 i *2 i *2 /o\ 

= a q Q + a xQ x + a y% + a zq z ■ ( 8 ) 

Our notation differs here from that often found in the literature (ao = t 2 ; a x = 
R 2 ] a y = R 2 ; a z = R 2 Z ). We have introduced the new notation since the fits to 
experimental data can give negative values for these parameters, especially if there are 
cross-terms. Throughout this paper we use the coordinate system in which the z axis 
is along the beam direction, the x axis is along the "out" [JTO | direction (the direction 



in the transverse plane along the momentum of the pair) and the y axis is along the 
"side" direction, which is perpendicular to the "out" and beam directions. Note that 
in the source rest frame the corresponding cross terms a x0 = a y0 = a z0 = 0, if the 
source velocity is defined as 0: 

(rd) - (ri)(t) , s 

*■= °- (9) 

We also assume that the source is azimuthally symmetric (in the source rest frame): 
CLxy = CLyz = 0, ci x = a y = a±, and there is no correlation of the u x-z" type: a xz = 0. 

As was mentioned above only three components of qi are independent, thus the 
experimental data can be fitted by the correlation function in the form: 

Q = A x q 2 x + A y q 2 y + A z q 2 z + A xz q x q z . (10) 

The correspondence of our parameters A to the often used R 2 is obvious, and we 
have adopted again a new notation for the reason given above. We assume that the 
source can move only longitudinally (we consider the more general case in the following 



section) and make use of the Lorentz invariance of the correlation function. After some 
algebraic manipulations one finds the relations between the parameters A, which are 
to be extracted from the data analyzed in the analysis frame defined by rapidity y, and 
the source parameters eij introduced above: 

A z = — -3— -[a z cosh 2 (y - y*) + a sinh 2 (y - y*)} (11) 

cosh (y-y) 

v± 2 

A x = a ± H -y— -[a cosh 2 (y - y*) + a z sinh 2 (y-y*)] (12) 

cosh {y-y) 

A y = aj_ (13) 
2V± 2 

A xz = — -3— r [a z cosh(?/ - y ) sinh(y - y*) 

cosh (y - y) 

+a cosh(y - y*) sinh(y - y*)] . (14) 
In the case when the analysis frame coincides with the source rest frame (y = y*) these 



equations (|lT|--[14|) give the well known equalities: 



A z = a z + a V?, (15) 

~ 2 

A x = a± + a V± , (16) 

Ay = Oy, (17) 

A xz = 2a V ± V z . (18) 

When the analysis frame coincides with the LCMS (y = y) we have: 

A z = a z cosh 2 (y - y*) + a smh 2 (y - y*), (19) 

A xz = 2V ± smh(y -y*)[a z - a ), (20) 

A x - A y = V ± 2 [a cosh 2 {y-y*)+a z smh 2 {y-y*)], (21) 



and it follows: 



A z - - - 2 y = a z -a , (22) 



A x A,, 
Vi 

A xz /V ± 



— = tanh(2(y-y)). (23) 
A z + (A x - A y )/V ± 



The last equation, in principle, permits to evaluate the velocity of the source with 
respect to the pair velocity. 

The equations presented above can be used as a check of the validity of the func- 
tional form (Gaussian) used for the correlation function. One needs to evaluate Ai as 
a function of analysis frame rapidity y and compare the observed dependence of the 
parameters on the rapidity with the dependence given by these expressions. One of 
the simplest tests of this type would be to look for the independence of the parameter 
A y (often called the "side radius squared") on the analysis frame rapidity y. 

In the discussion above it is assumed that the source rest frame exists, that is that 
all pions entering the detector are emitted by a source with fixed rapidity. In general 
this is not true. It is rather likely that the effective source should be described as a 
superposition of sources moving with different rapidities. In this case, approximating 
the correlation function by a single (multidimensional) Gaussian will not be successful; 
it could result, for example, in the dependence of A y , the effective source "side" size, 
(and/or the coherence parameter A) on the analysis frame rapidity. This arises from 
the fact that the radii A4 (other than A y ) depend on the analysis frame rapidity y 
(see Eq.^4|). If the source parameter a y is different for the different sources, a fit to 
a single Gaussian results in an apparent ^/-dependence for the A y . Note that in the 
saddle point approximation [|| the effective pion source is exactly in the (Gaussian) 
form considered above. In this case the question about the existence of the source rest 
frame is equivalent to a question about the validity of the saddle point approximation. 

3 Directed flow effects 

Here we discuss how directed flow affects the terms in the correlation function related 
to the "side" direction. In particular we argue that in the presence of directed flow the 
"side-long" and "side-out" cross terms of the correlation function become non zero. 
The experimental measurements of these terms (especially as a function of centrality 
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of the collision) would be of great interest. 

We start with the simple case when the correlation function can be written in the 
"source rest frame" in a pure Gaussian form (|8|-|T0|) without any cross terms except 
the trivial ones related to the fact that particles are on the mass shell. We assume also 
that directed flow exhibits itself simply as a movement of the source in the transverse 
plane. Under this assumption one can calculate the correlation function in the analysis 
frame using Lorentz transformations. To do this we start with the expression for the 
correlation function in the source rest frame(H-fTD|) , then perform a Lorentz shift in the 
direction of the transverse flow (with transverse rapidity y±), and then do a Lorentz 
shift along the z axis to the frame moving with rapidity y. In the final step we use the 
mass shell constraint, Eq. ([[]). We denote by ip the angle between the direction of the 
source velocity (the reaction plane) and the ("out" direction). 

This gives the expressions: 



A * = TYT- \ cosh 2 (y-y*) + a sinh 2 (y - y*) cosh 2 y ± 

cosh (y - y) 

+a± sinh 2 (y — y*) sinh 2 y±), (24) 

. . 2 , / , i V± cosh(?/ — y*) sinhy_|_ 2 

A x = a j. sin ip + a_i_ (cos ip cosh y± H — ) + 

cosh(y - y) 

V±cosh(y-y*)coshy ± . 2 V± smh 2 (y-y*) 

+a ( T-rz r h cos V; sinh y ± ) + a z -j- — ,(25) 

cosh(y — y) cosh (y — y) 

A y = aj_ (cos 2 ip + sin 2 ip cosh 2 y±) + a sin 2 ip sinh 2 y±, (26) 

a n ■ i i, f i v. V ± cosh{y -y*)smhy ± 

A xv = 2aj_ sin^ cosh y± ( cos ip cosh y± H — r J 

cosh(y — y) 

. . . V± cosh(?/ — y*) coshyj_ 

+2a sin ip sinh y±_ ( — h cos ip smh y± ) 

cosh(y - y) 

-a ± sin(2^), (27) 

a o-i, sinh(y -y*) f V± cosh (y - y*) sinh y ± 

A xz = 2a± smh y ± — (cos ip cosh y ± H — ) 

cosh(y — y) cosh(?/ — y) 

smh(y-y*) V± cosh(y - y*) coshj/ ± 

+2a coshyx — -(cos^smhT/i. H — ) 

cosh(?/ — y) cosn(y — j/J 
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+2a z V ± sin% - y») roS ffi_ y *{ (28) 
cosh (y - y) 

A ( . ,sinh(2y ± )sinV>sinh(y-?/*) 

= (a± + ao) r-rr^ r • (29) 

smh(y — y) 

Let us consider, for example, the "side-long" cross term A yz . Averaging it directly 
over the reaction plane angle ip gives zero, but the expression one has to average is not 
A yz , but the correlation function exp(— A yz q y q z ). Then directed flow effectively results 
in nonzero cross terms in the correlation function. Note that analogous terms appear 
also in A x , A xy , A xz . 

If no selection of events is done with respect to the reaction plane orientation, 
averaging over ip gives: 

— j dipexp(-A yz q y q z ) = I (B yz q y q z ), (30) 



By. = (a ± + a ) &miH ^r Ky ,~" J , (31) 



where 

( sinh(2yj_) sinh(y — y*) 
sinh(y - y) 

and one could fit the data using the corresponding form of the correlation function. 

We should mention that strictly speaking the averaging over the flow angle should 
be performed separately in the numerator and denominator in the definition of the 
correlation function, Eq. (0), rather than averaging the correlation function itself. The 
above result is valid if the width of the two particle correlation related just to the 
movement of the source (this correlation has nothing to do with the quantum inter- 
ference of the identical particles) is much wider than the width of the correlation due 
to quantum interference. This problem/assumption is relevant not only to anisotropic 
flow effects. One should not forget that the HBT correlations are always "on top" of 
any long and/or short range correlations due to resonance decays, source movement, 
energy and momentum conservation, etc.. By fitting the correlation function in a form 
appropriate only for quantum interference correlation we implicitly assume that all 
other correlations have much greater widths. 
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Formulae (p4H29D can also be useful for purposes not directly related to transverse 
flow. For example, the case of ip — and tanh(?/_|_) = V± cosh(y — y*) would give the 
parameters of the correlation function as measured in the pion pair rest system. 



4 The ratio of the correlation functions 

We now suggest a method to study the functional form of the correlation function, 
exploiting its Lorentz transformation invariance. First we hypothesize a particular 
functional form of the correlation function. We then calculate the parameters in dif- 
ferent Lorentz systems and compare the observed dependence on the rapidity of the 
analysis frame with the expected dependence. In other words we choose the form and 
parameters of the correlation function and test for its Lorentz invariance. 
In practice we suggest the fitting of the ratio of the correlation functions: 

dNSL 

R(a x ,a y ,a z ,a ;X,v x ,v z ) = — (32) 

where (1) and (2) denote different event subsamples and/or different Lorentz systems. 
Note that in this method there is no need to produce mixed pairs; it also much less 
affected by the corrections of various types (such as the Coulomb correction). 

To illustrate the method we consider an example related to directed flow. One 
can assume that the source of pions in events with different reaction plane orientation 
differs only in orientation of directed flow velocity, but all other parameters describing 
the source are the same. Then one can perform a fit to extract the source transverse 
velocity. One should fit the ratio: 

. 1 + Xe- Q ' . . 

R{a x ,-) = 1 + Ae - Q > ( 33 ) 

where Q! and Q differ only by a Lorentz shift in the flow direction. For example, if one 
can select two event subsamples with the flow direction pointing along the +x direction 



and the — x direction, then Q' and Q would be defined by formulae (24-29) with cosip 
equal to +1 and —1, respectively. 
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One could argue that the source might look different from the +x and — x direc- 
tion [|7| , that is the parameters of the effective sources emitting pions in these directions 
are different. In this case it could be better to use event subsamples corresponding to 
flow directed along the +y and — y directions. In this case, because of symmetry, the 
only difference between Q' and Q is due to the direction of the movement of the source. 

Note that a simple test of whether the parameters of the source are the same 
for different event subsamples (selected, for example, according to the orientation of 
the pair momentum with respect to the reaction plane) could be the evaluation of 
the correlation function using invariant variables. The simplest case would be a one 
dimensional fit in Qi nv . If the parameters of the source for different event subsamples 
are the same, and the only difference between them is the direction of the motion of 
the source, the extracted R inv should be the same for all event subsamples, due to the 
Lorentz invariance of the correlation function. 

5 Likelihood function fitting 

The usual approach for the fit of the correlation function assumes the introduction of 
binning in some part or even over the entire phase space. This approach has a number 
of disadvantages. It is sensitive to the experimental acceptance and to cuts introduced 
during the analysis. Often some bins have a low number of hits, which introduces 
biases in the results. Sometimes the parameterized correlation function has only a few 
parameters, but the complicated functional form of the parameterization may demand 
binning in many dimensions. 

It is attractive in this case to use the maximum likelihood fitting technique, which is 
free from almost all of the problems mentioned above. However a direct application of 
this technique is difficult, because in this approach one needs to parameterize not only 
the correlation function itself (the goal of the fit), but of the distributions of "true" 
and "mixed" pairs. The latter is extremely difficult, since one must take into account 
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numerous nontrivial acceptance and analysis cuts. 

We propose a method based on the likelihood function, which allows one to avoid 
these problems. For simplicity we consider one dimensional case, but the generalization 
of the method for multiple dimensions is straightforward. In the one dimensional case 
the correlation function is written as: 

c{x) K ( — )/( ^^ } - (34) 

Let us denote the probability density for mixed pairs as p(x): 

1 dN m i xec i p{x^j (35) 

N m ixed dx 

Then the probability density for the true pair distribution can be written as 



1 dNf rue 



= G[C]C(x)p(x), (36) 



N true dx 

where G is the normalization factor, which is a functional of C(x), and therefore 
depends on the parameters which determine C. 

The log-likelihood function depends on the correlation function in the following 

way: 

lnL[C] = Yl {lnG[C]+\nC(x k )+\np(x k )}. (37) 

k={true pairs} 

The last term in the sum does not depend on C (x) and can be omitted. But one must 
evaluate the normalization factor G[C] for different sets of parameters. The problem is 
that for this calculation one needs to know the probability density p(x). The trick of our 
method is that the required integral can be evaluated using "mixed" pair distribution. 

G-^C] = f dxC(x)p(x) = f dxC(x) 1 dN ™ xed (38) 

J J iv mixed 

« E C( Xi ) (39) 

mixed i = { m i xe d pairs} 

For example, if we consider the Gaussian form of the correlation function 

C(x) = 1 + Aexp(-ax 2 ), (40) 
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the log-likelihood function can be written as 

In L(a, A) = £ ln(l + Aexp(-ax^)) (41) 

k={true pairs} 

-N true \n{j^— J2 (l + Aexp(-ax?))}. (42) 

mixed i = { m i X ed pairs} 

The method has many advantages. First, it gives a correct evaluation of the errors in 
the fit, because it automatically takes into account the fact that the pairs in the mixed 
event technique are not independent and in general the fluctuations are proportional 
to 5N mixed oc N^lf xed [|TTJ, in contrast to the fluctuations in true pairs distributions 

1/2 

vN true oc N t r Ue . Taking into account that the maximum number of mixed pairs which 
can be generated is proportional to the square of the number of events (N mixed oc 
^events 00 ^true)i ^ follows that the fluctuations in mixed pairs distributions never 
become much less than that in the real pairs distribution, regardless of how many 
of them are generated. Other advantages of the likelihood technique is that it is 
free from the problems associated with finite bin size, and it is very flexible in its 
applications (For example, it is easy to introduce analysis or acceptance cuts, and it is 
very straightforward to use different functional forms of the correlation function.) 

6 Conclusion 

We have presented new techniques for the analysis of like-particle pairs using the HBT 
method and for testing the form of the correlation function used in such studies. Start- 
ing from the Lorentz-invariant property of the correlation function, we derive relation- 
ship for the fitted parameters as a function of the rapidity of the analysis frame. This 
technique is particularly useful, but not limited to, studies of interacting systems in 
which anisotropic transverse flow is present. 

Our results may be summarized as follows: (i) The analysis of the same data in 
different Lorentz systems helps in the extraction of the source parameters and assessing 
the validity of the approximations used, for example the validity of Gaussian source 
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(saddle point approximation), (ii) Anisotropic transverse flow would result in the non 
zero "side" cross terms in the correlation function, (iii) The study of ratios of the 
correlation functions evaluated in different Lorentz frames (the ratios of "true" pair 
distributions) is an important new analysis tool, (iv) We have presented a new fitting 
technique based on the likelihood function which does not require binning and has 
many advantages in comparison to the conventional method. 
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